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Q ' We consider a particle which moves on the x axis and is subject to a constant 

00- 

force, such as gravity, plus a random force in the form of Gaussian white noise. We 

analyze the statistics of first arrival at point xi of a particle which starts at xq with 
O ' velocity vq. The probability that the particle has not yet arrived at xi after a time 

^ ! 

^ ' t, the mean time of first arrival, and the velocity distribution at first arrival are 

C^ , all considered. We also study the statistics of the first return of the particle to its 

-t— > . 

starting point. Finally, we point out that the extreme-value statistics of the particle 
and the first-passage statistics are closely related, and we derive the distribution of 



n^ • the maximum displacement m = maxt[x(t)]. 
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I. INTRODUCTION 

In this paper we consider a particle which moves on the x axis and is subject to both 
a constant force, such as gravity, and a random force in the form of Gaussian white noise. 
The Newtonian equation of motion is given by 

(.,(t)) = 0, {,,(()#')> = 2Ai(i -(') , (2) 

where g is a. constant. 

Simple stochastic processes such as this are of both mathematical and physical interest. 
The approximately random collision forces experienced by a particular particle in a many 
particle system are often modelled by Gaussian white noise. Langevin's equation [l|] for 
the motion of a Brownian particle in a constant force field corresponds to Eq. ([1]) with an 
additional viscous damping term of the form —A dx/dt on the right-hand side. On setting 
(7 = in Eq. ([1]) and regarding t as a Cartesian coordinate instead of time, one may interpret 
the path x{t) of the particle as a configuration of a semi-flexible polymer [2|. For several 
applications of the process ([1]) related to semi-fiexible polymers and driven granular matter, 

n n n n 

see [21, 13|, U, 15[ and references therein. 

In this paper we study first-passage properties [Gj of the process ([1]). More precisely, we 
analyze the statistics of the first arrival at point Xi of a particle which starts at xq with 
velocity vq. Due to translational invariance no generality is lost in choosing xi to be the 
origin, and since we consider both positive and negative g, no generality is lost in choosing 
Xo to be positive. Thus, "first passage" corresponds to the first exit of the particle from the 
positive X axis. If the initial velocity Vq is positive, the particle must return to its initial 
position at least once before exiting from the positive x axis. Thus, in the limit Xq \ with 
f > 0, the first-passage statistics reduces to the statistics of first return of the particle to its 
initial position. Clearly, first-passage statistics, as defined here, is the same as the statistics 
of absorption of a particle moving on the half line a; = with an absorbing boundary at 
x = 0. 

First-passage properties of the random acceleration process, corresponding to Eq. ([1]) 
3ut without the constant term q on the right hand side, are derived or reviewed in Refs. 
ifl fl fl a. In .he .e.ai„d. of tM= .c*„ we .„w .ow U..e ..1. ca„ be gen.al.ed .o 



include the constant force. In Section [TTl some statistical quantities of interest in connection 

our explicit results are presented. In Section 



with first passage are defined, and in Section [IT 
IIVI we show that the extreme- value statistics 4, llO, lUl, ll2, ll3| and first-passage statistics 
of the process ([T]) are closely related. This is then used in deriving the distribution of the 
maximum displacement m = max([x(t)] of a particle which begins at the origin with velocity 

Vq. 

With no loss of generality we replace the parameters g and A, introduced in Eqs. ([1]) 
and ([2]), by g — > 7 = ±1 and A = 1 throughout this paper, since this can be achieved by 



rescaling [1J| the variables x and t. For 7 = — 1 and 7 = 1, the constant force drives a 
particle on the positive x axis toward and away from the origin, respectively. 
Integrating the equation of motion ([T]) yields 

x{t) =xo + vot + l 7^2 + {t- t')r]{t')dt' , (3) 

Jo 

which, together with properties ([2]) of the random force, implies the moments 

(x) =xo + vot + \^e , ((x - {x)f) = 1^ . (4) 

Thus, the contribution of the random force on the right side of Eq. ([3]) has typical size t^/^. 
For large t the constant force is more important than the random force, but for small t the 
opposite is true. 

It is convenient to define P^{x,v;xo,vo;t)dxdv as the probability that the position and 
velocity of a particle, moving according to Eq. ([1]) with (? — > 7 = ±1 and A = 1, evolve from 
Xq, Vq to values between x and x + dx, v and v + dv in a time t without ever reaching x = 0. 
The probability distribution Pj satisfies the time-dependent Fokker-Planck equation [l| 

d d d d"^ 

h V h 7 

dt dx dv dv"^ 
with the initial condition 

P^(X, V] Xo, Vq] 0) = 5{X - Xo)5{v - Vq) (6) 

and the boundary condition 

P^(0,t;;xo,t;o;i) = 0, v>Q. (7) 

This boundary condition ensures that only trajectories which exit the positive x axis for the 
first time at time t contribute to P-y(0,t;; xo,fo; ^)- Trajectories which leave the positive x 



^ + ^^r: + ^37 ~ ^;3 ) ^7(^'^5a;o,t'o;t) = 0, (5) 



axis at an earlier time and return to the positive x axis are excluded. Equation ([7]) is also 
the appropriate boundary condition for motion on the half line x > with an absorbing 
boundary at x = 0. 

In the absence of the constant force, the corresponding probability distribution 
Po{x,v]Xo,vo;t) satisfies the same Fokker-Planck equation ([5]), initial condition ([6]), and 
boundary condition ([7]), except that the term '-fdP/dv in Eq. ([5]) is absent. This implies the 
relation 

Py{x, v; Xo, Vq; t) = exp [^jiv - Vq) - \t\ Pq[x, v; Xq, Vq; t) (8) 

between the distributions with and without the constant force, which is central to our work. 
In a classic paper on the first-passage properties of a randomly accelerated particle, McK- 
ean [7(] derived the exact propagator Po(0, —v; 0,fo; t) for t> > and Vq > 0, corresponding 
to a particle which leaves the origin with velocity vq and returns for the first time at time t 
with velocity —v and speed v. His result and Eq. (jHD imply 

P^{0,-v;0,vo;t) 

= ^ exp [-i 7(^ + Vo)-lt~{v'- Wo + vl)/t\ erf U^^ , (9) 

where erf(z) denotes the standard error function [ISl. Il6l|. 
The Laplace transform 

POO 

P^{x,v;xo,Vo;s)= / dte~'^^P^{x,v;xo,Vo;t). (10) 



plays a central role in our work. Substituting Eq. (^ on the right-hand side, using the 



integral representatio n |15l . Il6| erf(z) = 2n ^^"^z J^^ dy exp(— z^y^), and integrating over t 
with the help of Ref . 16| , we obtain 

P^{0,-v;0,Vo;s) 

2vr Jo 

X [{v^ - Wo + vl + ?,woy^)-^'^ + (4s + lY'^{v^ - wq + v^ + 3woyY^] . (11) 

We will also need McKean's result ^] for the Laplace transform, 

P^{0,-v;0,Vo;s) 

= ^ / df,^^^^KU^A^^TT)v)K,,{yqAFTT)vo), (12) 

TT Wfo Jo cosh(^7r/i) 



where K^{z) is a modified Bessel function [l5|, [l6|. Expressions (ITT]) and (fT2|) are particularly 



convenient for numerical and analytical calculations, respectively, and both expressions are 
used below. 

The exact solution Pq{x,v\xq,vq\ s) of the Fokker-Planck equation for random accelera- 
tion on the half line x > with boundary condition ([7]), is given in Ref. [2], where it is 
derived from more general results of Marshall and Watson [8]. All of our results for first 
passage from an arbitrary initial point Xq are based on this solution. Substituting it in Eq. 
(IH]) and setting x = 0, we obtain 

/■oo 

P^(0, -V- xo, v^- s) = e-^(-+''«)/2 / dFe-^^Vs+i/4,F(-t^)^s+i/4,F(^o) , (13) 

Jo 

where 

^P^^piv) = F-i/^Ai [F^/h + F~^'h) , (14) 

and Ai(2;) is the Airy function [l5|. Some important properties of the two set of basis 
functions V's,f(^) and 0s, ^('y) are discussed in Ref. |2[. For example, (ps^piv) vanishes 
identically for i; > 0, so that Eq. (iT3l) satisfies the boundary condition (171). 



II. STATISTICAL QUANTITIES OF INTEREST 

The " survival probability" or probability that a particle with initial position and velocity 
Xo and fo has not yet left the positive x axis after a time t is given by 



oo poo 



Q^{xo,vo;t) = I dv dx P^{x,v;xo,vo;t) . (16) 

-oo Jo 



According to Eqs. fl5l). fiTl). and flT6l) 

d f°^ 

— Q^(xo, vo; t) = - dv vP^{0, -v; xq, vq; t) . (17) 

Thus, we interpret 

vPy{0,—v;xo,vo;t)dvdt, (18) 

for Xq > 0, as the probability that the particle reaches the origin for the first time at a time 
between t and t + dt with speed between v and v + dv. 
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Several useful relations follow from this interpretation of the quantity (fTSl) . The survival 
probability defined by Eq. flT6|) . its limiting value for t —>■ oo, and its Laplace transform can 
be written in the form 

ft 



pT noo 

Q.y{xQ,vo;t) = 1 — / dt' / dvvPj{0,—v;xo,vo;t'), (19) 

Jo Jo 

/■oo 

Q^{xo,vo;oo) = 1 - / dvvP^{0,-v;xo,Vo;0) , (20) 

Jo 



Q-y{xo,vo;s) = - 



'0 

fOO 



/■oo 

1— / dvvP^{0,—v;xo,Vo;s) 
Jo 



(21) 



The mean time to exit the positive x axis for the first time is given by 

/o°° dt t /q°° rft; t;P-,(0, -v; xp, vp; t) 
jp dt J^ dvvPy{u,—v;xo,Vo;t) 

= -[1-Q^{xo,vo;oo)y^ Urn— dvvP^{0,-v;xo,Vo] s) . (23) 

Since a particle which begins at the origin with a negative velocity immediately moves 
onto the negative x avis, Q^ and T^ satisfy the boundary conditions 

Q^{0,vo;t) = 0, t;o<0, t>0, (24) 

T^{0,vo) = 0, vo<0. (25) 

Finally, the probability that the speed of the particle on exiting from the positive x axis 
for the first time is between v and v + dv is given by G^{v; Xo,Vo) dv, where 

^ / s vJ^dtP^{0,-v;xo,vo;t) vP^{0,-v;xo,vo;0) ._. 

Jq dvv Jq dtP^{0,-v;xo,Vo]t) 1 - Q^[xo,vo;oo) 

and the normalization f^ dvG^{v; xq, vq) = 1 has been imposed. 



III. RESULTS 

A. Limit Q^{xq,vq;oo) of the Survival Probability 

The survival probability Q^{xo, Vq; t) introduced in the preceding section can, in principle, 
be evaluated for arbitrary t from Eqs. ([9]), flT3l) . and flT9|) . However, this involves integrating 
over t' and v and, for xq > 0, inverting a Laplace transform to go from s to t, most of which 
must be performed numerically. In this section we consider the limiting value Q^{xo,vo; C)o) 



or probability that in an infinite time the particle never leaves the positive x axis, which can 
be obtained analytically. 

In the absence of a constant force, Qo{xo,vo,t) decays as t~^^^ in the long-time limit 
[21, lZI, Isj. Thus, Qo{xo^vo, oo) = 0, which means that the particle exits from the positive x 
axis in an infinite time with probability 1. As shown in Eq. ([3]), the constant force adds an 
extra term ^ 7^^ to the displacement x{t) of the randomly accelerated particle. In the case 
7 = —1 of a constant force toward the origin, the particle leaves the positive x axis sooner 
than without the constant force, so 



Q-i{xo,vo;oo) = 



(27) 



In the case 7 = 1 of a constant force pushing the particle away from the origin, the 
corresponding probability Qi{xq,vq,oo) does not vanish and is expected to increase as Xq 
and fo increase. Recalling Eq. I^^ . substituting Eqs. ( iT2l) and ( iT3l) into Eq. (l20l) . and 
proceeding as described in the Appendix, we obtain 



I , t;o < , 
Qi(0,t;o;oo) = <( . r- — 

erf I w 1 170 



Wo >0 



(28) 



Qi(2;o,t;o;oo) = 1 



-yo/2 



27r 



y" dF F~'l^ exp (- Y^ - Fx^ Ai {F^'\, + \ F'^/^) 



(29) 



As in Eqs. ([9]) and (IT4|) , erff^:) denotes the error function, and Ai(z) is the Airy function, 
both defined as in Ref. 15|]. We found it useful to make the change of variables F = u~^ in 
evaluating the integrals in Eqs. ( l29l) . ( l34l) . and ( l50l) numerically with Mathematica. 
Equations fl28l) and fl29|) imply the asymptotic behavior 



Qi(0,fo;oo) 



'' '« \ 1/2 

, Vo\0 

71 / 



1/2 



(30) 



(5i(xo,0;oo) 






~3i>o/2 



xo \0 



f — > CXD 



^TT^Xo 



1/4 



-(2x0/3)1/2 



(31) 



Xo ^ CXD 



Qi(a;o,t;o;oo) 



erf 



33/2 xo 



d 



Qi{x,vo;oo) 



/2^|t;op/2 



dx 

i>o/2-|i>oP/9a;o 



in'^Xo 



1/4 



-t)0-{2a:o/3)i''2 



J x=0 

vo<0 

Xq 



Xo , Vo>0 
, a;o \ , 



Xo 



(32) 



In Fig. 1, Qi{xo, Vq] t), as given by Eqs. (!28|) and (!29|) . is plotted as a function of Vq for 
several values of Xq. Note that (5i(xo,fo;i) increases monotonically with increasing xq and 
Vq, as expected. 

If the random force is switched off, the particle trajectory becomes x(t) = xo + vot+ \'~it^i 
which never reaches the origin for 7 = 1 and Vq > —^Ix^. Thus, each of the smooth curves 
in Fig. 1 is replaced by the unit step function Qi(xo,fo; C)o) = Q{yQ + ^2xq). The short 
vertical lines in Fig. 1 indicate the value of Vq at which Q\{xq, Vq\ 00) is discontinuous. 



B. Mean First-Passage Time T^(xq,vq) 

The t~^/^ decay of the survival probability Qq{xq, fo, t) of a randomly accelerated particle 
[2I, lZI, Isj is so slow that the mean time of its first exit from the positive x axis, given by 
To(xo,fo) = ^^ dtt\—dQQ{xQ,VQ,t^ld'i\, is infinite. However, in the case 7 = —1 of a 
constant force toward the origin in addition to the random force, the corresponding mean 
time T_i(a;o,'yo) is finite. Recalling Eq. ([25]), substituting Eqs. ([12]), ([13]), ([28]), and ([29]) 
into Eq. (1231) . and proceeding as described in the Appendix, we obtain 



, t;o < 



T-i(0,t;o) = <^ W— e 



2^0 _-„o/2 



vr 



2 



2 — erfc 



(33) 



_ e^''°/2erfc(v^) , t'o > 
o 



1 As" /"°° 

T_i{xo,vo) = -J- dt t-^/2 (g^,^ + 2vQt + t^) exp 



3 f xp + vpt _ 1 ,/, 



t3/2 



-t^ 



3^0/2 
~2^ 



j^ dF F~"' exp (- Y^ - i^xo) Ai (FV3,;, + 1 F-^/^) 



X 



^-^^f-G^-bb 



(34) 



where erfc(2;) = 1 — erf (2;) is the complementary error function 15|. 



Equations (133!) and (1M|) imply the asymptotic behavior 



18^0^^/^ 



T.,iO,v,)^{ \ n J ' ^0^0' (35) 

2f + I , t;o -^ oo , 

T_,(.o,0)-<| W ''''F(t)"°''' "°^°' (36) 

The leading terms 2fo and (2xo)^^^ in Eqs. ( 135|) and (136!) for xq = 0, I'd -^ oo and for 
fo = 0, Xo — > cxD are easy to understand. According to the discussion following Eq. (jlj), 
the constant force is more important than the random force for large t, implying x{t) ^ 
xo + vot — ^ t^, which vanishes at T_i(xo, vq) ^ Vq + a/wq + 2xo. 

For small xq and vq the particle tends to reach the origin quickly, so T^(xo,fo) is small. 
For short times the random force is more important than the constant force (see discussion 
following Eq. (jl])) and primarily responsible for the asymptotic behavior T-j-i(0,fo) ~ v^^'^ 
and T_i(xo, 0) ~ x^^^ in Eqs. (!35|) . (l36l) . and (138|) . We note that these same power laws for 
small Xq and vq appear in the mean first exit time [171, [l8[ of a randomly accelerated particle 
with initial position and velocity xq, fo from the finite interval < x < L. 

In Fig. 2, T_i(xo,fo), as given by Eqs. ( l33l) and ( IMl) . is plotted as a function of vq 
for several values of Xq. As expected, T_i(xo,fo;i) increases monotonically as Xq and Vq 
increase. 

As discussed just above Eq. ( 128|) . for 7 = 1 the probability Qi(xo, "^0; 00) that the particle 
never leaves the positive x axis is nonzero, in general, and it increases, as in Eqs. (l30ll and 
( 13T|) . as fo and xq increase. However, the mean first exit time for those trajectories which 
do leave the positive x axis, defined by Eq. (|23|) . is finite. From Eqs. (IT2|) . (|23|) . (I25l) . and 
( 128|) . we obtain 



/ 



which has the asymptotic behavior 

Ti(0,wo;oo 



, t;o < 

2 f 

3 W TT 3/ g^f^/^./s^ 



2 _^ f /6^ 2 \ e-^-°/^ (37) 



o \ 1/2 



37r 



2.0-^^ +3, .0-00 



(38) 
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For arbitrary xq and Vq, Ti(xo,fo) follows from substituting Eqs. (fT3|) and (129|) into Eq. 
f l23|) . This leads to a lengthy expression, containing multiple integrals, which we were unable 
to simplify and omit here. 

The results ( 133|) and (1371) for T_i(0,fo) and Ti(0,fo) are compared in Fig. 2. 

C. Distribution G^{v;xo,vo) of the Particle Speed at First Passage 

For arbitrary initial position xq and initial velocity vo, the distribution G^{v]Xo,vo) of 
the particle speed on exiting from the positive x axis for the first time is determined by Eqs. 
flT3|) . fl26|) . and fl29|) . Here we restrict our attention to the case of a particle which begins at 
the origin with vq > 0. The distribution G^{v; 0,fo) of its speed on returning to the origin 
for the first time and leaving the positive x axis can be readily evaluated by substituting Eqs. 
(ITTl) and (l28ll into Eq. ( l26l) and integrating over the variable y numerically. The results, for 
several values of vq, are shown in Fig. 3. 

Each of the curves in Fig. 3 has a single peak. As vo increases, the peak shifts to larger 
values of v, as expected, and becomes broader. The peak position and width correspond 
roughly to the mean speed (v) at first return and the root-mean-square deviation a = 
(f - (f )^) ) , where 



K). 



dvv^G^{v] 0,fo) . 



(39) 



Using Eqs. flT2l) . fl26l) . and fl28l) and the approach of the Appendix, we have calculated 
the first two moments of the speed at first return analytically, obtaining 



{v)_ 



Vo + i: + e 



X 



2v^ 

71 



1/2 



-yo/2 



vo + - ] e'^^/^erfc 






V, 



+ 2t;o + 2(t;)_, 



/ 



-3^o/2 



V 



erfc 



(^')i = 4 - 2vq + 



/ 

2i^ 
37r 



1/2 



(^0 - 5) 



- - e^^^erfc (v^) 



:(^)i + l 



(40) 
(41) 

(42) 
(43) 



For large vo the average speed at first return and the root-mean-square deviation from the 
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average have the asymptotic forms 



(v)^ 



a 



and 



(v). 



(^1 



Vo + 



-Vq 

3 



1/2 



(44) 
(45) 



Us- 27r)i;o 



.1/2 



(46) 
(47) 



which are quahtatively consistent with the evolution of the curves in Fig. 3 as Vq increases. 

For large Vq the constant force is more important than the random force, and x(t) ~ 

Xo+Vot+^'jt'^. Thus, for 7 = — 1, the particle returns to its starting point with approximately 

the same speed it had initially. This is consistent with the asymptotic behavior in Eq. (jH 



IV. EXTREME- VALUE STATISTICS 



Consider a particle which begins at xq = with velocity vq and moves according to Eq. ([T]) 
with f? — > 7 = ±1. At some time in the interval < t < 00 the particle attains a maximum 
displacement m = maxt[x{t)]. For large t, x(t) ^ |7^^, as follows from the discussion below 
Eq. (jl]). Thus, in the case 7 = 1 of a constant force in the positive direction, m = 00. In this 
section we consider the less trivial question of the maximum displacement m for 7 = — 1, 
and we derive the corresponding distribution V_-\(, m, Vn) ■ D istributions such as this play a 



central role in the field of extreme- value statistics 
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ll| . The extreme- value statistics of 



a generalized Gaussian process that includes random acceleration as a special case is studied 
inRefs. Il2,ll3j. 



To derive the distribution V-i{'m,Vo), we begin by writing 

d 

V-i{m,Vo) = -— J'_il(m,t;o) , 
am 



(48) 



where JF_i(m, fo) is the probability that, for a constant force in the negative direction, the 
displacement x{t) of a particle which begins at the origin with velocity Vq never exceeds m 
in the time interval < t < 00. For m < 0, jF_i(m, t>o)=0, since the initial displacement 
xq = already exceeds m. For m > 0, jF_i(m, vq) is the same as the probability that, for a 
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constant force in the positive direction, a particle with initial position m and initial velocity 
—vq never reaches the origin. This follows from the invariance of the probability under the 
coordinate transformation x ^ m — x. Since this latter probability is precisely the survival 
probability Qi{m, —vq; oo) considered in Sections HTl and UTTl 

T-i{m, vq) = 9{m) Qi{m, -vq; oo) , (49) 

where 6{m) is the standard step function. 

Making use of Eqs. (HHll and (H9l) and the expressions for Qi{0, vq; oo) and Qi{xq, vq; oo) 
in Eqs. ([28D and ([29]), we obtain 



V-i{m, Vq) = 6{-Vo) erf ( ^| \vo\ j 6{m) 

„Do/2 /-oo / 1 \ 

+e{m) -y= / dF F-^/^ exp ( - Fm J Ai {-F^/^o + \ F'^/^) (50) 

for the extreme- value distribution. The distribution vanishes for m < and is normalized so 
that j_ dmV-i{m, vq) = 1, as follows from Eqs. ( HHI) and ( 149|) and the boundary condition 
Qi(oo, —Vq] oo) = 1. The first term on the right side of Eq. fl50|) has its origin in the non-zero 



probability erf ( ^ /| \vo\ ) (see Section [III A|) that a particle which begins at the origin with 



f < never returns to the origin, in which case the maximum displacement m equals the 
initial value xq = 0. 

The extreme- value distribution V^iijn, vq) is plotted as a function of m for several positive 
and negative values of Vq in Figs. 4a and 4b, respectively. In the absence of the random 
force, X = fot — |t^, which implies P_i(m, t>o) = 9{—vo)6{m) + 9{vo)S{m — ^VQ). The random 
force broadens the delta functions, as seen in the figure. 

For positive Vq, the peak in Fig. 4a shifts to larger values of m and becomes broader 
as Vo increases, as expected. The mean value and the root-mean-square deviation vary as 
(m) ^ ^Vq + Vq and a ~ (f'^o) for large positive vq. For large vq the constant force is 
more important than the random force, and the leading term in (m) equals the maximum 
displacement |wq of a particle subject only to the constant force. 

In the results for iiq < in Fig. 4b, the vertical line at m = represents the term pro- 
portional to 5(r7i) in Eq. fl50|) . The most probable value of m, which maximizes V-i{m, Vq), 
is zero for all negative vq- The mean value of m is positive and, for vo negative and large in 
magnitude, (m) ^ |(2|t;o|/37r)i/2e-3l''ol/2^ and (m^) ^ f (2|t;o|V37r)i/2g-3|t,o|/2_ 
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V. CONCLUDING REMARKS 

This completes our study of the first-passage and extreme- value statistics of the process 
([T]). In closing we note that in a mathematical tour de force, Marshall and Watson (8| derived 
the Laplace transform Pg,\{x, v;xo,vo]s) of the solution to the Klein-Kramers equation with 
the absorbing boundary condition ([7]). The Klein-Kramers equation is the Fokker-Planck 

equation for the process |l[ 

d X dx 

_ + A- =« + „(*), (51) 

which, unlike Eq. ([T]), includes viscous damping and plays a central role in the theory of 
Brownian motion. In principle, all of the first-passage and extreme- value properties we have 
considered follow, for this more general process, from the expression of Marshall and Watson 
for Pg,\{x, f ; Xo, fo; s), which, however, involves an infinite double sum over special functions 
and is difficult to work with. 
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APPENDIX A: CALCULATIONAL DETAILS 

The quantities Q^{0,vo;oo) and T^(0,fo), defined in Eqs. (!20l) and (l23l) . can both be 
expressed in terms of the integral 

oo 

dv vP^{0, —V] 0, vq] s) 



^-7»;o/2 roo sinh(7r/x) , , , r- , ,, , 

Jo cosh(i7r/i) 



tt'^Vq Jq ■ ■ coslil^l 



oo 



F^(/i, s) = / dve-''^/''K,,{,/iAs + l) v) , (A2) 



where we have used the expression for P^(0, —v; 0,fn; s) in Eq. (IT2|) . Evaluating F^(/i, s) 



with the help of the integral representation [15|, ll6 1 

/■oo 

Ki^{y)= dtcos(/it)e-''™^^* (A3) 

Jo 
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and substituting the result in Eq. (lAip . we obtain 



e 



-71)0/2 



U,2\-l/2 



rff t;P^(0, -w;xo,t;o;s) = (4s + 1- h"^) 

f^ sinh r/iarccos(i7(4s + 1)^^/2)1 , ,.- , , , 

x/ ci/x/i ^^ A^ZL^ ^ ^ir,^(v/(47TI)t;o). (A4) 

Jo cosii (gTr/zj 

First we consider the quantity Q^(0, Vq; cxd), defined in Eq. fl20|) . For 7 = — 1 and s = 0, 
the arccos in Eq. flA4p equals |7r. From Eqs. (!20l) and (\AAh and the relation [16| 

/■OO 1 

/ rf/i/i smh{bfi)Ki^{vo) = -TTVosmbe-^^^"'^ , (A5) 

Jo ^ 

we obtain Q_i(0,t'o; 00) = for vq > 0, in agreement with Eq. fl27|) . 

For 7 = 1 and s = 0, the arccos in Eq. (]A4p equals |7r, and expression (128|) for 
Qi(0,fo; C)o), with fo > 0, follows from Eqs. ( !20l) and ( ]A4p and the relation 

/oo AT / /o \ 

d/i/i tanh (ivr/i) i^i^(t;o) = ^ 7rwoe''°/^erfc I V ^ 1 ' (^6) 

which we derived with the help of the integral representation ( JA3I) . 

We now turn to Q-y{xo,Vo; 00) for xq 7^ 0. The result Q-i{xo,vq; 00) = was established 
in the paragraph containing Eq. flTTI) . Expression fl2^ for (5i(a;o) ^o! 00) follows from Eqs. 
flT3|) - flT5|) and fl20|) . The lengthy derivation will not be given here, but it is easy to see, with 
the help of Ref. [2i], that the result (129|) satisfies the appropriate Fokker-Planck equation 
{vodx^ + -yd^^ + dl^)Q^{xQ, Vq; cxd) = and to check, by numerical integration, that Eqs. fl28|) 
and ( !29l) agree for a;o = 0. 

The results ([33D and (1571) for T^(0, fo) may be derived by substituting Eqs. (1251) and (IA4p 
into Eq. (125]) and using Eqs. (lASp . (1A6P and some analogous integrals over /x which can be 
evaluated with the help of the integral representation (\A3\i . Expression (IMI) for T_i(xo,fo) 
follows from substituting Eqs. (IT3|) - (IT5l) and (|27 [ ) into Eq. (I23l) . but the derivation is long 
and will not be given here. Making use of Ref. 2|, we have confirmed that the result (134|) 
satisfies the appropriate Fokker-Planck equation {vodx^ +'^9^^ +d'^^)Ty{xo,vo; 00) = —1 and 
checked by numerical integration that Eqs. ( I33l) and ( l34l) agree for xq = 0. 

The moments (f") , defined by Eqs. ( 1260 and (15^ . can be expressed in terms in terms 

of the integral 

f°° ~ p-i^ol-i / r5 \ " r 1/9 

I dvv-+'P,{Q,-v-x,,v,-s) = ^——[-2^\ |(4. + 1-172)-^/^ 

[^ sinh [u arccos (i7(4s + l)-^/2^1 ^ ^ 

Jo cosh (3vr/ij 
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This relation is the same as Eq. (lA4p except for the extra factor f " introduced by applying 
{—2d/d'~^Y to the quantity F^{fi,s) in Eq. flA2l) . The steps leading from Eq. flA7l) to the 
results for (t>)^ and (f^)^ in Eqs. ( H0l) - (H3l) are very similar to the steps, described above, 
from Eq. ( ]A4[) to the final expressions for Q^{0,vo; oo) and T^{0,vo). 
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FIG. 1: Probability Qi{xq, vq; oo), given by Eqs. (128]) and (f29|) . that a particle subject to a random 
force plus a constant force pushing it away from the origin never leaves the positive x axis. If the 
random force is switched off, Qi{xo, vq; oo) becomes a unit step function, as discussed in the second 
paragraph below Eq. ([5T|) . The short vertical lines in the figure indicate the values of vq at which 
the step function jumps from to 1 for xq = 0, 0.1, 1. 
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FIG. 2: (a) Mean time T_i(xo, vq) to exit the positive x axis for the first time for a particle subject 
to a random force plus a constant force toward the origin, given by Eqs. (j33p and (j34p . (b) Mean 
first exit times r_i(0, fo) and Ti(0, wq), given in Eqs. ([33]) and ([37j) for constant forces toward and 
away from the origin, respectively. 
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FIG. 3: Distribution G^{v;0,vo), given by Eqs. (fTTj) . ([25]) . and ([25]) . of the particle speed t; at 
first return, for vq = 0.2, 0.4, 0.8, and 1.6. The results in (a) and (b) are for constant forces 
directed toward and away from the origin, respectively. As vq increases, the peak becomes lower 
and broader and moves to the right. 
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FIG. 4: Distribution V-i{m,VQ), given by Eq. (j50p . of the maximum displacement m attained by 
a particle which begins at the origin with velocity vq and moves according to Eq. ([T]) with g = —1. 
The results in (a) are for vq = 0.4, 0.6, 0.8, and 1.0. As vq increases, the peak becomes lower and 
broader and moves to the right. The curves in (b) correspond, from top to bottom, to vq = -0.4, 
-0.6, -0.8, and -1. The vertical line at m = represents the term in P-i(m,VQ) proportional to 
5{m). 



